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Abstract  A  model  for  oscillating  free 
surface  jet  flow  of  a  fluid  from  an 
elliptical  orifice,  together  with 
experimental  measurements,  can  be 
exploited  to  characterize  the  elonga- 
tional  viscosity  of  non-Newtonian 
inelastic  fluids.  The  oscillating  jet 
flow  is  predominantly  elongational, 
with  a  small  strain  that  oscillates 
rapidly  between  large  and  zero  strain 
rates.  We  find  that  to  reproduce  the 
experimentally  observed  steady 
oscillating  jet  flow  in  model  simula¬ 
tions,  the  assumed  form  of  the  non- 
Newtonian  viscosity  as  a  function  of 
strain  rate  must  have  zero  gradient, 
i.e.,  be  Newtonian,  at  zero  strain  rate 
(a  behavior  exhibited,  in  general,  by 
real  inelastic  fluids).  We  demonstrate 
that  the  Cross,  Carreau,  Prandtl- 


Eyring,  and  Powell-Eyring  forms, 
although  they  have  finite  viscosity  at 
zero  strain  rate,  have  either  nonzero 
or  even  unbounded  gradient  at  zero, 
and  hence  are  unable  to  model  oscil¬ 
lating  jet  behavior.  We  propose  a  new 
non-Newtonian  viscous  form  which 
has  all  of  the  desirable  features  of 
existing  forms  (high  and  low  strain 
rate  plateaus,  with  adjustable  loca¬ 
tion  and  steepness  of  the  transition) 
and  the  additional  feature  of  Newto¬ 
nian  behavior  at  low  strain  rates. 
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surface  jet  •  Strain-rate  thinning 
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Introduction 

When  a  fluid  exits  an  elliptical  orifice  it  establishes  a  jet 
with  a  free  surface  that  is  steady  in  time  and  oscillates  in 
space  (Fig.  1).  Surface  tension  is  the  restoring  force  and 
viscosity  is  the  damping;  in  broad  strokes,  the  wave¬ 
length  of  the  oscillating  free  surface  is  indicative  of  the 
surface  tension  of  the  fluid  in  the  ambient  atmosphere, 
and  the  decay  of  the  oscillation  indicates  the  fluid’s 
viscosity.  Lord  Rayleigh  (1879)  developed  a  relation 
between  surface  tension  and  the  wavelength  of  oscilla¬ 
tion  for  an  inviscid  jet  in  the  absence  of  gravity  which, 
when  combined  with  experimental  measurements  of 
oscillating  jets,  determines  the  surface  tension  of  the 
fluid.  Since  then  oscillating  free  surface  jets  have  been 
used  to  measure  surface  tension;  the  advantage  of  this 


technique  is  that,  because  the  surface  of  the  jet  is  newly 
created,  the  technique  measures  dynamic  surface  ten¬ 
sion,  i.e.,  tension  as  a  rapidly  changing  function  of  time 
for  newly  formed  surfaces.  Pedersen  (1907),  Bohr  (1909), 
Hansen  et  al.  (1958),  Defay  and  Hommelen  (1958),  and 
Thomas  and  Potter  (1975)  each  present  improvements  to 
Rayleigh’s  relation.  These  algebraic  equations  for  sur¬ 
face  tension  in  terms  of  the  measured  wavelength  and 
amplitude  of  the  free  surface  oscillation  neglect  gravity, 
are  based  on  a  small  departure  of  the  oscillating  profile 
from  a  circular  cylinder,  and  model  the  fluid  as 
Newtonian  viscous.  In  previous  work  (Bechtel  1989; 
Bechtel  et  al.  1995,  1998)  we  derive  models  for  a  finite 
departure  from  equilibrium  which  account  for  gravity; 
these  models  take  the  form  of  nonlinear  integro- 
differential  equations  for  the  free  surface  profile.  In 
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Fig.  1  The  oscillating  jet  phenomenon;  from  Lord  Rayleigh  (1890) 

Bechtel  et  al.  (1995)  we  employ  the  inverse  formulations 
of  these  models,  in  conjunction  with  experiments 
consistent  with  the  model  assumptions,  to  determine 
Newtonian  viscosity  as  well  as  surface  tension. 

The  model  equations  employed  in  all  of  the  inverse 
techniques,  except  Bechtel  et  al.  (1995),  assume  a 
constant  surface  tension  and  constant  viscosity  within 
the  wavelength  over  which  the  measurements  are  taken. 
Bechtel  et  al.  (1995)  generalize  the  model  to  allow  for 
variable  surface  tension  and  a  non-Newtonian  viscous 
fluid  (although  the  non-Newtonian  capability  is  not 
exploited  in  that  paper).  In  Bechtel  et  al.  (1998)  we 
investigate  the  signature  in  Newtonian  oscillating  jet 
behavior  of  surface  tension  which  changes  within  a 
wavelength  and  on  submillisecond  timescales.  Here 
we  investigate  the  effect  of  non-Newtonian ,  strain- 
rate-dependent  viscosity  on  oscillating  jet  behavior. 

In  the  investigation  we  have  discovered  that  existing 
forms  of  strain-rate  dependent  viscosity  (Cross,  Carreau, 
Prandtl-Eyring,  and  Powell-Eyring)  all  exhibit  anoma¬ 
lous  non-zero  gradients  at  low  strain  rates,  so  that  they 
are  unable  to  capture  the  Newtonian  plateau  exhibited  in 
general  in  real  inelastic  fluids  at  low  strain  rates;  Analysis 
reveals  that  the  viscosity  gradient  at  zero  strain  rate  is 
unbounded  for  the  Cross  form,  and  finite  but  nonzero  for 
the  Carreau,  Prandtl-Eyring,  and  Powell-Eyring  forms. 

In  physical  experiments,  slowly  varying  (i.e.,  slender) 
oscillating  jet  flows  of  strain-rate-thinning  non-Newto¬ 
nian  viscous  fluids  are  readily  observed,  so  the  3-D  free 
surface  boundary  value  problem  should  have  solutions 
which  are  slowly  varying.  The  kinematics  of  oscillating 
jet  flow  are  such  that  if  the  effect  of  gravity  is  negligible 
(as  is  often  the  case),  the  strain  rate  is  zero  at  the  extremes 
of  each  free-surface  oscillation.  When  the  Cross,  Car¬ 
reau,  Prandtl-Eyring,  and  Powell-Eyring  forms  are 
employed  in  the  slender  jet  equations,  they  each  fail  at 
these  locations  of  zero  strain  rate  because  of  their  non¬ 
zero  gradients  there,  erroneously  predicting  the  develop¬ 
ment  of  small  length  scales  and  subsequent  jet  breakup. 
The  discovery  that  the  slowly  varying  asymptotics  can 
produce  no  consistent  solution  indicates  a  failure  of  the 
constitutive  model  for  the  fluid;  the  mathematical  model 
has  incorrect  assumptions,  and  the  only  empirical  ad  hoc 
model  feature  is  the  viscosity  form. 


Our  observations  have  an  important  ramification  in 
the  modeling  of  non-Newtonian  viscous  fluids  in  general. 
The  oscillating  jet  flow  serves  as  a  test  case.  If  a  proposed 
viscosity  form  cannot  model  the  oscillating  jet  flow,  it  is 
suspect  to  be  applied  to  any  flow  that  passes  through  low 
strain  rates.  An  acceptable  non-Newtonian  viscous 
constitutive  form  must  yield  a  solution  for  the  oscillating 
jet  problem  which  is  slowly  varying  and  stable  when  the 
fluid  being  modeled  is  observed  to  undergo  a  stable, 
slowly  varying  oscillating  jet  flow;  to  do  this  the  form 
must  possess,  in  addition  to  a  bounded  zero  strain  rate 
viscosity,  a  sufficiently  small  gradient  at  zero  strain  rate, 
i.e.,  it  should  be  effectively  Newtonian  at  low  strain  rates. 

Motivated  by  this  inadequacy  of  existing  character¬ 
izations  of  non-Newtonian  inelastic  fluids,  we  present  a 
new  viscosity  form,  called  the  Zhou  form,  that  has  a  zero 
strain  rate  viscosity  gradient  which  is  identically  zero, 
i.e.,  it  is  Newtonian  at  low  strain  rates.  We  find  that  the 
Zhou  form,  unlike  the  other  forms,  leads  to  slender 
oscillating  jet  solutions. 


The  governing  equations  for  a  steady  viscous 
oscillating  jet 


In  Bechtel  et  al.  (1995)  time-dependent  governing  equa¬ 
tions  for  an  oscillating  slender  jet  of  a  viscous  fluid  are 
derived  allowing  for  variable  viscosity  and  surface 
tension.  We  recall  only  essential  features  of  this  deriva¬ 
tion,  which  reveal  why  the  various  strain-rate  dependent 
viscosity  forms  are  inconsistent  with  slender  oscillating 
jet  behavior. 

The  fluid  is  assumed  incompressible,  and  flows  in  the 
direction  of  gravity.  We  define  Cartesian  coordinates 
such  that  the  e3  direction  coincides  with  the  centerline  of 
the  jet,  and  assume  the  jet  cross  section  is  elliptical  for 
several  oscillations,  so  that  the  free  surface  is 
2  2 

_ N _ a - — _ i  =  o  .  (1) 

$7  (*3,1)  ®l(x3,t) 


The  fluid  we  model  is  viscous  (not  viscoelastic),  with  a 
strain-rate-dependent  viscosity.  Since  the  oscillating  jet 
flow  is  predominantly  elongational,  the  appropriate 
viscosity  is  elongational  viscosity,  which  we  label  ty  if 
the  fluid  is  Newtonian,  ;;  is  constant  and  three  times  the 
shear  viscosity.  Here  the  viscosity  is  variable,  an  explicit 
function  of  the  second  invariant1  of  D: 


i/(Hd)D,  Hz)  - 


tr(D2)  -  (trD)2 


(2) 


where  T  is  the  determinate  part  of  the  Cauchy  stress 
tensor  T  —  pi  and  D  is  the  symmetric  part  of  the  velocity 
gradient.  For  an  incompressible  fluid  (trD  =  0),  1\D  is 


'The  definition  (2)2  is  not  universal,  e.g.,  Bird  and  Armstrong 
(1987)  define  their  second  invariant  II  =  4/;-(D2),  so  that  II  =  8IID. 
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non-negative.  The  customary  definition  of  strain  rate 
magnitude  y  is 

y  =  2 y/Uo  .  (3) 

In  order  to  deduce  the  dominant  physical  effects  and 
exploit  the  slenderness  (or,  equivalently,  slow  variation) 
of  the  jet,  the  problem  is  non-dimensionalized  through 
adoption  of  characteristic  scales  r0,  z0,  f0,  v0  for 
transverse  length,  axial  length,  force,  and  axial  velocity, 
respectively.  The  dimensionless  forms  x,  y,  z,  x,  0i,  02  of 
the  coordinates  jc1}  x2,  x3,  t  and  evolving  semi-axes  <3>! ,  ®2 
of  the  elliptical  jet  cross  section  are  introduced  through 

.  zo 

Xi  =  r0x,  x2  =  r0y,  x3  =  z0z,  t  =  —x, 

(4) 

T>i(x3,t)  =  roxjf^z,  t);  ®2(x3,t)  =  ro02(z,l)  . 

The  small  parameter  in  the  slender  jet  theory  is 
the  slenderness  ratio  e  =  ^  <C  1.  For  an  oscillating  jet  e 
is  roughly  the  ratio  of  the  mean  radius  of  the  jet  cross 
section  to  the  wavelength  of  the  oscillation,  typically  0.1. 
The  velocity  field  within  the  oscillating  jet  is 

v  =  r0{[£xCi(z,r)  +0(e3)]e!  +  [eyC2(z,x)  +  0(e3)]e2 
+(v(z,x)  +  0(e 2))e3}  ,  (5) 


where  C\,  £2,  v  are  dimensionless  0(1)  functions;  note 
that  the  flow  is  elongational.  For  this  flow  the  second 
invariant  (2)2  is 


Ho  —  ^2  |Ci  +  C2  +  ^  >  (6) 

where  “,z”  denotes  differentiation  with  respect  to  the 
dimensionless  axial  coordinate  z.  To  leading  order  in  the 
slenderness  ratio,  II/j,  and  hence  the  viscosity,  are 
functions  only  of  the  axial  coordinate  and  time. 

To  the  fixed  observer  of  an  oscillating  jet  experiment, 
the  free  surface  profile  does  not  change;  given  the  focus 
of  this  paper  on  variable  viscosity,  we  restrict  ourselves  to 
constant  surface  tension.  Therefore  we  recall  the  form  of 
the  governing  non-dimensional  equations  for  the  oscil¬ 
lating  jet  which  is  steady  with  constant  surface  tension: 


y  1  ,Z  y  I  i  , 

4l  —  T — j  <=2  —  T — 1  V(P\(P2  —  1) 

01  02 

BP  -  2Zct f»j 02Ci  +^41Kc  +  s2^p-  (Z_-CliZ  +  ZCliZZ) 


0102 


W 

2\ 


(vCi  ,z  +  Ci)) 


BP-2Z^^2C2  + 


thKs+£2  Ml 

W  4 


(z,zC2jZ  +  zc2jZZ) 


e2  (0102  ,zKs  +  0201, z  Kc)  +  £2(2  0i02Zzrz 
+  20102Zu^z  -BP)  T  ~  0i  02  =  rr,z0!02 


(V) 


where 


D  _  Z5/o  1  _  ffZ0  1  _  9Z 0  ry,^  _  f/(Z)Z0 

2  4  )  tjz  23’zr  2  ^  lzl 

(r  oosin  r  ir 


7  ,  ,  h,z(z)z0 


0 


3provo 


3pr2v0  ’ 

LJJ  (P  -  Pa)&C\dX2, 


P  = 


nfo 


r2n 


Kc  =  - 
n 


k  cos2  8d6 


(8) 


0102 


r»27r 


cos2  8 


71  0°  (02  sin2  9  +  03  cos2  0) 


3/2 


d0  +  0(e2), 


r*27r 


**  =  - 
71 


k  sin2  Odd 


0102 


r*27r 


sin2  8 


o  sin"  8  +  0;  cos2  0) 


3/2 


d0  +  O(e2)  ; 


p,  a,  g,  pa  are  the  fluid  density,  surface  tension, 
acceleration  of  gravity,  and  specified  ambient  pressure, 
assumed  constant,  and  k  is  the  dimensionless  curvature 
of  the  surface. 

Equations  (7)!_3  are  the  leading  order  equations 
which  follow  from  the  kinematic  free  surface  boundary 
condition  and  incompressibility.  These  equations  are 
independent  of  the  constitutive  nature  of  the  fluid  and 
hold  for  all  regimes  of  jet  behavior.  The  constant  of 
integration  in  the  flow  rule  (7)3  presupposes  that  the 
characteristic  length  r0  is  the  geometric  mean  of  the 
principal  radii  of  the  elliptical  cross  section  at  z  =  0  and 
the  characteristic  velocity  v0  is  the  axial  velocity  at  z  =  0. 

Equations  ( 7)4-6  correspond  to  projections  of  con¬ 
servation  of  momentum  in  the  x,  y,  and  z  directions, 
respectively,  and  are  special  to  the  non-Newtonian  viscous 
constitutive  assumption  (Eqs.  2).  We  have  retained  the 
leading  order  terms  within  each  of  the  viscosity  (Z  and 
Zz),  surface  tension  (IF),  constraint  pressure  (B),  gravity 
( F ),  and  inertia  ( 1 )  terms  (if  surface  tension  is  non¬ 
constant,  IFZ  terms  appear  in  Eqs.  (l)^;  see  Bechtel 
et  al.  1998).  For  a  Newtonian  fluid  the  dimensionless 
parameter  Z  z  is  identically  zero.  The  magnitude  of  the  full 
set  of  dimensionless  parameters  (8)1  _5  relative  to  the 
slenderness  ratio  s  indicates  the  relative  balances  of  the 
competing  physical  effects  in  the  free  jet  problem.  In  a 
particular  experiment  some  subset  of  these  effects  will  be 
dominant,  and  thus  survive  in  the  leading  order  equations. 
The  jet  oscillates  if  inertia,  surface  tension,  and  the 
constraint  pressure  are  leading  order  in  the  transverse 
momentum  equations  (7)4  5;  this  condition  translates  to 
jf  B «  0(1).  We  set  B=  1  by  fixing  the  force 

scale  /o  =  ^3^,  after  which  three  independent  scales 
remain:  r0,  z0  ^typically  10ro  for  slender  jets),  and  v0,  which 
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depend  on  the  experimental  conditions  and,  together  with 
p,c 7,  and  g,  completely  specify  the  oscillating  jet  behavior. 

The  consequences  of  variable  viscosity  in  the  asymptotic 
balance  of  Eqs.  (7)^  can  be  deduced  by  studying  three 
tenns: 

L0  =  IZcfrMU  LX  =E-^1\Z^  +  H^ZZ\, 

L2  =  2e2 1 Z-v-  +  Zv-z  -  P\  .  (9) 


L0  is  the  leading  order  viscosity  term  and  L\  the  next 
order  correction  in  the  transverse  momentum  equations 
(7)4,5,  and  L2  is  the  leading  order  viscosity  term  in  the 
axial  projection,  (7)6.  LQ  must  be  0(1)  for  oscillating 
viscous  jets  (if  LQ  is  0(e2)  or  less  the  oscillating  jet  is 
inviscid,  with  no  damping  to  leading  order).  The 
equations  (7)4_6  simplify  to  slender  jet  equations  only  if 
the  correction  terms  L\  and  L2  are  0(s2);  otherwise  these 
terms  compete  in  the  leading  order  equations  with  the 
formally  leading-order  L0  term,  thereby  disordering  the 
perturbation  expansion. 

To  anticipate,  in  the  analysis  immediately  following 

we  show  that  for  the  power  law,  Cross,  Carreau,  Prandtl- 

Eyring,  and  Powell-Eyring  forms  giving  viscosity  //  as  a 

function  of  second  invariant  11D,  the  viscosity  gradient 

.d'1  is  nonzero  at  IIn  =  0,  i.e.,  at  zero  strain  rate.  This 
d(iiD) 

feature  is  critical  in  the  oscillating  jet  characterization 
precisely  because  in  that  flow  the  strain  rate  oscillates 
about  zero:  we  discover  in  the  computations  of  the  next 
section  that,  since  (0)  ^  0,  the  dimensionless 

viscosity  gradient  Z  z  becomes  too  large  at  the  axial 
locations  of  small  strain  rates  for  these  standard  forms, 
causing  L\  to  be  huge.  This  implies  that  either  the  jet  flow 
cannot  remain  slender  or  the  empirical  viscosity  forms 
are  problematic.  Since  there  is  no  experimental  evidence 
of  an  inability  to  perform  oscillating  slender  jet 
experiments,  we  conclude  these  empirical  viscosity  forms 
have  unphysical  features  for  small  strain  rates. 

We  define  the  non-dimensional  strain  rate  magnitude 

y  by 


5  =  ^05 

Vo 


z0  ■ 
x— y 
2v0 


(10) 


The  dimensionless  form  of  the  viscosity  //(IId)  is  then 


Z  =  Z(~f)  •  (11) 

To  leading  order  in  the  oscillating  jet, 

~y  =  \J\(t2x  +  tl  +  vi)  .  (12) 


Power  law  fluid 


In  the  power  law  model  (Bird  et  al.  1987)  viscosity  is 
described  by  the  expression2 

2In  Bird  et  al.  (1987),  r\  =  kJ^  ,  which  is  equivalent  to  Eq.  (14), 
noting  II  =  8IId  from  the  previous  footnote. 


h  =  *(4IIz#  ,  (13) 

where  the  consistency  index  k  and  power  law  exponent  / 
are  positive  constants.  The  fluid  is  strain-rate  thinning 
for  l  <  1,  Newtonian  for  1=1,  and  strain-rate  thicken¬ 
ing  for  /  >  1. 

For  the  power  law  form  the  viscosity  gradient  is 

535)  =  *('-■  (i4) 


The  dimensionless  forms  of  the  viscosity  (Eq.  13)  and 
viscosity  gradient  (Eq.  14)  are 


z  =  zr?-‘,z'^J±flzP?- > 

ay  4 

with  the  two  dimensionless  constants 


(15) 


l,  ZP  = 


kzn 


(2vq 


1- 1 


,  2  i  I  •  (16) 

ipr^vo  \  z0 

For  a  strain-rate  thinning  fluid,  as  the  strain  rate  goes  to 
zero, 


>l{ 0)  oo, 


dry 


d(IIi 


■(0)  — >  —oo,  Z(0)  — »  oo,  Z(0)  — >  —oo  . 

(17) 


Cross  model 


The  predictions  for  a  strain-rate  thinning  fluid  of  infinite 
viscosity  when  the  strain  rate  is  zero,  and  zero  viscosity 
when  the  strain  rate  is  large,  are  shortcomings  of  the 
power  law  model.  Cross  (1965)3  proposed  a  model  with 
two  additional  constants  ;;0  and  r\ ^  that  corrects  these 
shortcomings: 

>i  =  n  oc+  I0-*1"  .  (18) 

1  +  (4*2IId)t 


If  /  <  1  and  t] o  >  >700  >  0,  the  viscosity  is  strain-rate 
thinning;  i]Q  and  are  the  limiting  viscosities  at  zero 
and  high  strain  rates,  respectively.  If  /  >  1  and 
ho  >  'loo  >  0,  the  fluid  is  strain-rate  thickening,  with 
>7o  and  r]^  reversing  roles  to  be  the  limiting  viscosities  at 
high  and  zero  strain  rates,  respectively.  The  Cross  form 
reduces  to  constant  viscosity  either  by  setting  1=  1 
(resulting  in  constant  viscosity  2fl+2oo)  or  setting  ?/0  =  t ]oa. 

The  viscosity  gradient  as  a  function  of  strain  rate  is 


d>/ 

d(HD) 


k'-1  ( 

ho  -  hoo) 

2'  (Hz# 

1  +  (4^11^)^ 

(19) 


3Cross  gives  >/  =  <7oo  +  jyyyj) j”>  in  which  D  =  2 flip.  This  is 
equivalent  to  Eq.  (18)  if  t  =  k,  m=  1-/.  In  his  original  paper  Cross 
fixes  m  =  I,  i.e.,  /  =  1. 
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The  dimensionless  forms  of  the  viscosity  (Eq.  18)  and 
viscosity  gradient  (Eq.  19)  are 


Z  —  Zoo  + 


ZO  -  Zoo 

l  +zcyl-r 


=  dZ(f)  _  {l  -  \){Z,  -  Zoo) 

df  C2y+i(\+Zc?~1)2  ’ 


(20) 


where  there  are  now  four  dimensionless  constants. 


/  Z 

1  5  -^OC 


zo 

o  hoo 

3provo 


Zo  = 


zp 


»7():  ZC 


2Au0\  1  ; 

.  zo  / 

(21) 


For  a  strain-rate  thinning  fluid,  as  the  strain  rate  goes  to 
zero  the  viscosity  is  bounded,  but  the  viscosity  gradient  is 
unbounded , 


'?(0)  =  '7o> 


d  i] 


d(IIz 


-(0)— >— oo,  Z(0)  =  Z0,  Z'(0)— >oo  ;  (22) 


the  Cross  fluid  is  strongly  non-Newtonian  at  low  strain 
rates. 


Carreau  model 


Another  remedy  for  the  power  law’s  unbounded 
viscosity  at  zero  strain-rate  is  the  Carreau  (1972) 
viscosity  model4,  which  has  been  shown  to  have  enough 
flexibility  to  fit  several  experimental  viscosity  curves 
(Bird  et  al.  1987;  Macosko  1993): 

'7  =  '7oo  +  (%-'7oo)(l+4A2IIz))¥  ,  (23) 

where  A  is  a  time  constant.  If  /  <  1  and  )/0  >  i]0 c  >  0 
this  form  models  strain-rate  thinning,  whereas  1  <  1  and 
0  <  ;/q  <  i] oo  models  strain-rate  thickening.  In  both 
cases,  j/o  and  17^  are  the  limiting  viscosities  at  zero  and 
high  strain  rates,  respectively.  The  form  reduces  to 
Newtonian  either  by  setting  I  —  1  (resulting  in  constant 
viscosity  r/0)  or  setting  »/0  =  rj ^ 

The  viscosity  gradient  as  a  function  of  strain  rate  is 

^=2A2(l-l)(rlo-rloo)(l+4A2lIDf  ,  (24) 


and  the  dimensionless  viscosity  and  viscosity  gradient  are 


z  —  Zoo  +  (Zo  —  Zoo ) 
Z'=Z2cy(l-  1)^ 


1  +  (Zcrj) 

z 


1  +  ( ZcyI) 


(25) 


where  the  four  dimensionless  constants  are  /,  the  Z0  and 
Zoo  defined  in  Eq.  (2 1  )o  3,  and 


4The  Carreau-Yasuda  model  is  given  in  Bird  et  al.  (1987)  in  the 
form  =  [1  +  (Ay)0]~°~  where  y  =  This  model  with  the 

parameter0 a  =  2,  referred  to  as  the  Carreau  model,  is  equivalent  to 
Eq.  (23). 


_  2Av0 

^ CY  ~  - 

z0 


(26) 


For  this  form,  both  the  viscosity  and  viscosity  gradient 
are  finite  at  zero  strain  rate  for  strain-rate  thinning 
fluids,  but  the  gradient  is  nonzero : 


l(0)  =  ho,  (°)  =  2 a2(1  -  l)(f/0  -  ri^, 

Z(0)=Z0,  Z'(0)  =z2CY(i-  i)(z°~Zoo:>  . 


(27) 


Again,  this  model  predicts  non-Newtonian  behavior  at 
low  strain  rates.  The  only  way  the  Carreau  form  can 
model  Newtonian  behavior  at  low  strain  rates  is  if  it 
models  the  fluid  as  Newtonian  at  all  strain  rates,  since 
the  viscosity  gradient  at  zero  strain  rate  is  zero  only  if 
/=  1  or  t] 0  =  rioo. 


Prandtl-Eyring  model 


One  of  the  first  viscosity  forms  obtained  by  a  molecular 
theory  was  the  Eyring  equation  (Brodkey  1967).  A 
modification  of  this  form,  the  Prandtl-Eyring  model 
(Skelland  1967)  has  been  used  to  describe  the  rheology 
of  strain-rate  thinning  fluids: 


>7(11  D)  = 


Br,  osinh-'^E) 
v/n ~d 


(28) 


where  B  is  a  constant.  The  limit  when  the  strain  rate  is 
zero  is  t] 0.  Like  the  power  law,  this  form  can  only  predict 
inviscid  behavior  at  high  strain  rates.  The  viscosity 
gradient  as  a  function  of  strain  rate  is 


d(nZj)  2  ,ifv/TTf 

which  is  finite  but  nonzero  at  zero  strain  rate  (as  with  the 
Carreau  form): 


_^_(0)  =  — *>- 

d(IIi>r  6^2 


(30) 


Powell-Eyring  model 


The  Powell-Eyring  model  (Skelland  1967)  adds  a 
constant  to  the  Prandtl-Eyring  form: 


>i(Ud)  =  ■ 


B(ho  -  hoo)  sinh 


v/IC 


(31) 


resulting  in  a  finite  limit  at  high  strain  rates,  but  the 
viscosity  gradient  is  unaltered  from  the  Prandtl-Eyring 
model,  with  the  same  zero  strain  rate  value  given  in 
Eq.  (30). 
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Zhou  model 

We  now  propose  a  form  which  has  a  finite  value  of 
viscosity  and  a  zero  viscosity  gradient  at  zero  strain 
rate,  representing  for  the  first  time  a  non-Newtonian 
viscous  fluid  model  which  is  truly  Newtonian  at  low  strain 
rates : 


Of  all  the  strain-rate  thinning  viscosity  forms  consid¬ 
ered,  we  will  find  that  because  of  its  zero  gradient  at  zero 
strain  rate  this  form,  and  only  this  form,  can  model  the 
oscillating  jet  flow. 

Example  characterization  of  a  conceptual  fluid 


'?0  +  *1  o 


'/o  -  n  c 


0} 


■/(Hr 


/(77z>)  =  (  Yf  )  exP{  B 


+  /(Hi 


(32) 


Specifying  A  >  0,  B  >  0,  l  >  0,  t/0  >  >  0  models  a 

strain-rate  thinning  fluid.  The  constants  rj0  and  are 
the  limiting  viscosities  at  zero  and  high  strain  rates, 
respectively,  A  is  the  strain  rate  at  the  center  of  the 
transition  from  the  low  strain  rate  viscosity  plateau  r]0  to 
the  high  strain  rate  plateau  t /<*,,  and  a  controls  the 
steepness  of  this  drop.  The  constant  B  is  proportional  to 
the  range  of  strain  rates  over  which  the  fluid  is  strictly 
Newtonian ,  i.e.,  the  viscosity  is  strictly  constant. 

The  viscosity  gradient  for  this  model  is 

-**-=(„  -n)  f{UD) 

d(IIzj)  ^°°  ?o) 


m>)=- 


IId 

A 


[1  +/(Hd)] 

B  ' 
+  Tk 


2  ’ 


exp<  B 


(33) 

In  dimensionless  form  the  viscosity  and  viscosity 
gradient  are 


Z  = 


1  -  ( ZHff^v\ 

^0  ^ oo  I 

[s 

M^] 

} 

2  l  +  (Z//y2)"exp<j 

T 

z'  =  (za 


exp  1 5  1  -  (ZHj2y  | 


1  +  (Z//")2)"exp|fi 


(34) 

where  there  are  now  five  dimensionless  constants:  a,  B, 
the  Zq  and  defined  in  Eqs.  (21)2,3,  and 


Zh  — 


Ac2 


(35) 


For  this  form  the  viscosity  is  finite  and  the  viscosity 
gradient  is  zero  at  zero  strain  rate  for  strain-rate 
thinning  fluids: 

'7(0)  =  ^  mL)  (0)  =  °’  z(0)  =  z°’  z(0)  =  0  ■  (36) 


To  demonstrate  the  differences  between  the  above  six 
non-Newtonian  viscosity  forms,  we  fit  all  of  the  forms  to 
the  same  set  of  elongational  viscosity  vs  strain  rate  data 
points,  the  circles  O  in  Figs  2,  3,  and  4.  For  complex 
fluids  elongational  response  cannot  in  general  be 
inferred  from  easy-to-perform  shear  measurements 
(Moore  and  Pearson  1975;  Baid  and  Metzner  1977; 
Chai  and  Yeow  1988;  Ramanan  et  al.  1997),  and  it  is  a 
significant  challenge  to  produce  viscosity  vs  strain  rate 
data  points  directly  from  elongational  measurements  of 
an  inelastic  fluid.  We  comment  that  the  inverse  formu¬ 
lation  of  the  mathematical  model  presented  here, 
together  with  experimental  measurements  of  the  oscil¬ 
lating  jet  free  surface  profile,  can  be  used  to  produce 
data  sets  like  the  circles  in  Figs.  2,  3,  and  4,  but  that 
remains  for  future  work.  For  the  demonstration  pur¬ 
poses  of  this  paper  we  have  created  a  data  set  which 
exhibits  strain-rate  thinning.  Explicitly,  we  create  a 
hypothetical  fluid  with  a  low  strain-rate  Newtonian 
elongational  viscosity  plateau  of  0.15  Pa  s  for  strain- 
rates  1  x  10“5  s-1,  1  x  10-4  s“\  1  x  10-3  s_1,  and 

1  x  10~2  s-1,  and  a  high  strain-rate  Newtonian  plateau 
of  0.006  Pa  s  for  strain  rates  2  x  102  s-1,  4  x  102  s-1, 
7  x  102  s_1,  and  1  x  103  s  -1.  The  transition  data  points 
at  the  strain  rates  0.02  s_1,  0.04  s-1,  0.07  s-1,  0.1  s-1. 


strain  rate  (s  ’) 


Fig.  2  Plot  of  log  of  viscosity  vs  log  of  strain  rate  y  for  strain-rate 
thinning  viscosity  forms  fit  to  the  experimental  data  (O):  power  law 
( — ),  Cross  ( — ),  Carreau  (■••);  coefficients  in  the  forms  are  given  in 
Table  1.  The  power  law  form  predicts  unbounded  viscosity  at  low 
strain  rates  and  a  viscosity  which  approaches  zero  at  high  strain  rates 
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Fig.  3  Plot  of  log  of  viscosity  vs  log  of  strain  rate  y  for  strain-rate 
thinning  viscosity  forms  fit  to  the  experimental  data  (O):  Prandtl- 
Eyring  Powell-Eyring  (■■•),  Zhou  ( — );  coefficients  in  the  forms  are 

given  in  Table  1 .  The  Prandtl-Eyring  form  predicts  a  viscosity  which 
approaches  zero  at  high  strain  rates 

0.2  s'1,  0.4  s'1,  0.7  s'1,  1  s'1,  2  s'1,  4,  s"1  7  s'1,  10  s'1, 
20  s'1,  40  s'1,  70  s'1,  and  100  s'1  are  generated  from  the 
Carreau  form  (Eq.  23)  with  ^  =  0.006  Pa  s,  i] o  = 
0.15  Pa  s,  2  =  7.57  s,  and  /= 0.068. 

We  fit  each  viscosity  form  to  the  data  as  follows:  the 
power  law  coefficients  are  obtained  using  a  least  squares 
fit  of  only  the  data  between  10“2  s'1  and  102  s'1,  the 
coefficients  in  the  Prandtl-Eyring  form  are  obtained  using 
a  least  squares  fit  of  only  the  data  between  10~5  s'1  and 
102  s'1,  and  the  coefficients  in  the  other  four  forms  are 
obtained  using  a  least  squares  fit  of  all  data  (even  though 
we  have  used  the  Carreau  form  to  generate  the  transition 
data,  we  fit  it  to  all  of  the  data,  since  as  we  will  see  the 
Carreau  form  is  incapable  of  capturing  the  Newtonian 
low  strain-rate  plateau).  As  is  customary,  the  fits  are  done 
on  log  > /.  The  resulting  characterizations  are  in  Table  1. 

Figures  2  and  3  give  the  customary  representation  of 
the  strain  rate  dependent  viscosity,  namely  a  plot  of  log 
of  viscosity  vs  log  of  the  strain  rate  y  =  The 

power  law  and  Prandtl-Eyring  forms  are  unable  to 
capture  the  nonzero  viscosity  plateau  at  high  strain 
rates.  All  of  the  forms  except  for  the  power  law  appear  in 
these  log-log  plots  to  have  a  Newtonian  plateau  at  low 


strain  rates.  This  appearance  is  deceptive:  non-log  plots 
of  viscosity  vs  strain  rate  and  viscosity  gradient  vs  strain 
rate  at  low  strain  rates  (Figs.  4  and  5,  respectively)  and 
the  left  column  of  Table  1  reveal  that  all  but  our  Zhou 
form  have  negative  viscosity  gradient,  rather  than  zero 
gradient,  throughout  the  low  strain  rate  “plateau.”  As 
can  be  seen  from  Figs.  4  and  5,  only  our  form  is  capable 
of  modeling  a  fluid  with  constant  viscosity  over  a  range 
of  low  strain  rates,  the  feature  that  proves  to  be  critical 
to  model  oscillating  jets.  Note  also  from  Fig.  3  that  our 
Zhou  form  is  capable  of  fitting  the  transition  behavior 
predicted  by  the  Carreau  form,  more  so  than  the  power 
law.  Cross,  Prandtl-Eyring,  and  Powell-Eyring  forms. 

We  now  apply  the  six  constitutive  forms  to  model 
oscillatory  jet  behavior. 


Viscous  oscillating  jets 


We  assume  that  the  material  properties  and  process 
conditions  are  such  that  W  and  Z  are  C(l)  and  Z'  and 
are  at  most  0(e2),  so  that  the  leading  order  equations 
from  Eq.  (7)  are,  in  order: 

y  4>l,Z  y  4*2, Z  i  j  , 

Cl  —  —7 —  5  4=2  —  —j—1  0102  - 
01  02 

P-2Z^+^Kc  =  ^(^  +  tt),  (37) 

1 

p-2ZC2  +  -as  =  ^(C2,z  +  C5),  v=i  . 

In  this  regime,  pressure,  surface  tension,  inertia,  and 
viscosity  balance  in  the  transverse  oscillation,  and 
gravity,  surface  tension,  viscosity,  and  viscosity  gradi¬ 
ents  are  too  weak  to  accelerate  the  jet  in  the  axial 
direction  in  this  order  of  approximation.  The  system  of 
Equations  (37)  reduces  algebraically  to  a  single  integro- 
differential  equation  for  the  free  surface  profile  c/q : 

(l  +  0*1 )  0 1  ,Zz  ~  ^01, z  —Kc)  +  16  Z0j0lz  =  0  , 

(38) 


where 

Ks-Kc 


0 ^  r2n  cos2  6  —  sin2  6  ^ 

71  Jo  (cos2  9  +  0^  sin2  6)  ~  ^ 


Table  1  Coefficients  in  the 
viscosity  forms  characterizing 
the  conceptual  fluid 


Model 

Coefficients 

Power  law 

k  =  0.0282  Pa  s048,  /  =  0.48 

-00 

Cross 

i/oo  =  0.15  Pa  s,  );0  =  0.006  Pa  s,  k 

=  6.18  s,  /  =  .11 

-oo 

Carreau 

i/oo  =  0.15  Pa  s,  )?(>  =  0.006  Pa  s,  X 

=  7.57  s,  /  =  0.068 

-15.4  Pa  s3 

Prandtl-Eyring 

r\0  =  0.15  Pa  s,  B  =  0.0301  s_1. 

-27.7  Pa  s3 

Powell-Eyring 

))o  =  0.15  Pa  s,  B  =  0.0231  s'1,  jy,*. 

=  0.006  Pa  s 

-46.8  Pa  s3 

Zhou 

rfa o  =  0.006  Pa  s,  t ;0  =  0.15  Pa  s,  A 
B  =  0.25,  a  =  1.9345 

=  1.4102  x  10'2  s'2. 

0 
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To  leading  order  the  non-dimensional  strain  rate 
magnitude  (Eq.  12)  is 


Note  that  in  this  regime  the  strain  rate  is  zero  at  all 
locations  where  —  0,  i.e.,  at  the  maxima  and  minima 
of  the  free  surface  oscillation  (recall  (j>i (z)  is  the  evolving 
semi-axis  of  the  elliptical  cross  section).  Hence,  in  this 
regime  of  oscillating  jet  behavior  the  strain  rate  repeat¬ 
edly  passes  through  zero,  twice  in  each  wavelength. 


When  the  power  law  expression  ( 1 5),  for  dimensionless 
viscosity  Z  is  inserted  into  Eq.  (38)  we  get 


<t>\ 


(40) 


Power  law  fluid 


Fig.  4  Non-log  plot  of  viscosity  vs  second  invariant  IID  for  the  strain- 
rate  thinning  viscosity  models  of  Figs.  2  and  3.  Experimental  data 
(O).  Only  the  Zhou  form  models  the  Newtonian  (constant  viscosity) 
behavior  of  the  fluid  at  low  strain  rates:  the  Cross  form  has 
unbounded  slope  at  zero,  the  Carreau,  Prandtl-Eyring,  and  Powell- 
Eyring  forms  have  negative  slope  at  zero,  and  the  Zhou  fonn  has  zero 
slope  at  zero 


The  correction  terms  defined  in  Eq.  (9)2j3  reduce  to 


2  a  4 


u  = 


L2  =  er 


K , 
W 


2(J,z)Z'  +  Zjz. 

2  Zy-^iy^  +  y2) 


(39) 


recall  that  a  solution  of  the  slender  jet  (equation  38)  is 
valid  only  if  these  terms  remain  small,  i.e.,  0(f). 


Fig.  5  Viscosity  gradient  vs  second  invariant  I1D  for  the  strain-rate 
thinning  viscosity  models 


(l  +  +  ~w~(Ks  ~~  Kc ) 

+  16Z>0^/sign(01^)|^>lz|/=  0  . 


(41) 


The  solution  of  the  dimensionless  integro-differential 
equation  (41)  for  the  free  surface  profile  depends  on 
initial  conditions  t/>i(0)  and  </>i>r(0),  and  three  constants 
W,  ZP,  and  /: 


=j>l(W,ZP,l)  .  (42) 

Figure  6  displays  leading  order  solutions  for  strain-rate 
thinning  (/  =  0.5),  Newtonian  (/  =  1),  and  strain-rate 


dimensionless  axial  distance  z 


Fig.  6  Free  surface  profiles  4>\ (4)  for  inviscid  jets  with  rate-dependent 
viscosity  given  by  the  power  law:  <j>l(0)=l.5,  <f> i  _(0)  =  0,  W=l, 
ZP=  0.032,  with  /  =  0.5  (strain-rate  thinning,  — ),  1=1  (Newtonian, 
. ),  1=  1.5  (strain-rate  thickening,  •••) 
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dimensionless  axial  distance  z 


Fig.  7  Dimensionless  strain  rate  magnitude  y  for  the  solutions  in 
Fig.  6.  The  strain  rate  is  zero  at  all  locations  of  extrema  of  the  free 
surface  oscillation 

thickening  (/=  1.5)  fluids,  holding  W  and  ZP  fixed. 
Figure  7  displays  the  accompanying  evolution  of  the 
strain  rate  magnitude;  these  simulations  confirm  that  the 
strain  rate  repeatedly  passes  through  zero  at  the  locations 
of  minima  and  maxima  of  the  free  surface  oscillation. 

We  now  check  if  the  correction  terms  L\  and  L2 
remain  0(e2)  everywhere,  as  must  be  the  case  if  the 
solution  of  Eq.  (41)  is  to  be  consistent  with  the 
assumptions  inherent  in  Eq.  (41).  L2  behaves  properly. 
Elowever  L\  is  in  general  unbounded  at  locations  of  zero 
strain  rate  for  all  0  <  /  <  1  (all  strain-rate  thinning 
fluids)  and  1  <  /  <  3  (a  range  of  strain-rate  thickening 
fluids).  From  Eq.  (16),  the  viscosity  gradient  Z'  goes  to 
infinity  for  the  power  law  fluid  as  strain  rate  y  goes  to 
zero  if  /  <  3;  therefore  the  product  (yz)~Z'  is  unbound¬ 
ed  at  zero  strain  rate  since  the  square  (yz)"  of  the  strain 
rate  gradient  is  not  zero  when  strain  rate  y  is  zero.  In 
addition,  y  is  not  zero  when  y  is  zero,  so  the  product 
Zy,z  is  also  unbounded  at  zero  strain  rate,  due  to  the 
unboundedness  of  viscosity  Z.  Consequently,  Lx  is 
unbounded  at  the  maxima  and  minima  of  oscillations. 
This  is  illustrated  by  the  numerical  computation  of  Lx 
shown  in  Fig.  8,  where  we  overlay  L\  on  the  jet  profile 
for  the  strain-rate  thinning  fluid  of  Fig.  6;  similar 
behavior  occurs  in  the  strain-rate  thickening  solution 
of  Fig.  6.  This  diagnostic  indicates  a  failure  of  the 
asymptotic  approximation  induced  purely  by  the  behav¬ 
ior  of  the  power  law  form  near  zero  strain  rate;  The 
isolated  spikes  in  L\  invalidate  any  power  law  solution 
of  the  leading  order  problem  as  a  leading  order,  slowly 


Fig.  8  Leading  order  free  surface  profile  (f> ](z)  of  Fig.  6  with  l  =  0.5 

(strain-rate  thinning,  . )  and  correction  term  Z-i/a2  ( — ),  showing 

spikes  at  locations  where  (j>i~  =  0.  The  spikes  in  C/e2  have  amplitudes 
in  the  numerical  simulation  between  95  and  5.5  x  1016 

varying  approximation  to  the  3-D  behavior  by  disor¬ 
dering  the  presumed  asymptotic  expansion.  The  spikes 
falsely  indicate  the  onset  of  small  length  scales  and 
breakdown  of  the  asymptotic  approximation.  The 
analysis  of  the  3-D  instability  of  an  oscillating  jet 
described  by  this  viscosity  form  is  beyond  the  scope  of 
this  paper,  but  is  irrelevant:  physical  fluids  are  observed 
to  establish  oscillating  jet  flows  without  the  onset  of 
small  length  scales,  and  hence  the  power  law  cannot 
characterize  these  fluids  in  this  flow. 

Cross  model 

Although  the  Cross  model  was  constructed  to  avoid  the 
problems  of  the  power  law  model  at  zero  strain  rate,  we 
again  encounter  (Fig.  9)  artificial  large-amplitude  spikes 
in  the  next  order  corrections  at  locations  of  zero  strain 
rate.  This  is  due  to  the  unbounded  viscosity  gradient 
gjjjL  (2!  in  dimensionless  form)  at  zero  strain  rate  in  the 
Cross  model  (see  Eq.  224). 

Carreau  model 

The  numerical  simulation  shown  in  Fig.  10  for  an 
oscillating  jet  of  a  strain-rate  thinning  fluid  with 
viscosity  modeled  by  the  Carreau  form  again  indicates 
large  spikes  in  the  correction  terms.  While  the  Carreau 
model  yields  viscosity  and  viscosity  gradients  that  are 
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dimensionless  axial  distance  z 

Fig.  9  Leading  order  free  surface  profile  </>\(z)  ( . )  and  correction 

term  Lx/e2  ( — )  for  an  oscillating  jet  with  strain-rate  thinning  viscosity 
given  by  the  Cross  model:  W=  1,  Z0  =  1.0,  Zao  =  0.0171,  Zc=276, 
and  /=  0.5,  showing  spikes  in  Lx  at  locations  where  <^l  z  =  0.  The 
spikes  have  amplitudes  between  63.4  and  9.7  x  1012 

finite  at  zero  strain  rate  (see  Eqs.  27)  so  that  Lx  is 
bounded,  the  magnitude  of  Lx  is  still  large  enough  to 
disorder  the  asymptotic  expansion.  This  is  because  the 
dimensionless  value  of  Z'  at  zero  strain  rate  in  the 
Carreau  form  is  proportional  to  the  difference  between 
the  dimensionless  zero  strain  rate  viscosity  Z0  and  large 
strain  rate  viscosity  Z^,  and  typically  this  difference 
is  large.  For  example,  Z'(0)  =  25,900  in  the  simula¬ 
tion  displayed  in  Fig.  10  (which  models  a  hypo¬ 
thetical  experiment  with  the  material  properties 
p  =  1.0  g  cm-3,  a—  24.2  dyne  cm-1,  t/0  =  0.0396  Pa  s,  f/oo  = 
0.0057  Pa  s,  2=  2.09  s,  and  /=  0.5,  and  process 
conditions  of  ro  =  0.05  cm,  ro  =  0.5  cm,  and  vo  =  220 
cm  s-1).  From  Fig.  10  we  see  that  this  large  value  of 
Z'(0)  produces  significant  spikes  in  L\. 

We  comment,  without  exhibiting  the  graphs,  that  the 
Prandtl-Eyring  and  Powell-Eyring  forms  give  results 
similar  to  but  slightly  worse  than  the  Carreau  form,  with 
large  spikes  in  Lx  at  the  locations  of  zero  strain  rate.  This 
is  to  be  expected,  given  that  they  demand  a  comparable 
but  slightly  larger  magnitude  of  viscosity  gradient  at  low 
strain  rates  (see  Table  1  and  Figs.  4  and  5). 

Zhou  model 

From  the  analysis  above  it  is  clear  that  we  can  control 
the  product  (y  )  Z'  in  Lx  at  zero  strain  rate  only  with  a 
constitutive  form  in  which  the  viscosity  gradient  at  zero 


dimensionless  axial  distance  z 


Fig.  10  Leading  order  free  surface  profile  ( — )  and  correction 

term  Li/e2  ( — )  for  an  oscillating  jet  with  strain-rate  thinning  viscosity 
given  by  the  Carreau  model:  W=  1,  Z0  =  0.036,  Zao  =  0.00541, 
ZCY  =  1839,  and  1=0.5,  showing  spikes  in  Li/r  at  locations  where 
0l  z  =  O.  The  maximum  amplitude  of  the  spikes  in  the  numerical 
simulation  is  40 


dimensionless  axial  distance  z 

Fig.  11  Leading  order  free  surface  profile  (j) i(z)  ( . )  and  correc¬ 

tion  term  Lfis2  ( — )  for  an  oscillating  jet  with  strain-rate  thinning  vis¬ 
cosity  given  by  the  Zhou  model:  W=  1,  Z0  =  0.4545,  Zoc  =  0.018, 
ZH=1  x  10~8,  a  =  1.81,  and  B  =  0.05.  Note  the  uniformly  small 
magnitude  of  the  correction  term,  indicating  the  leading  order 
solution  is  a  valid  representation  of  the  3-D  solution 


383 


strain  rate  is  zero.  The  Zhou  form  (Eq.  34)  is  construct¬ 
ed  expressly  to  have  this  property.  The  simulation  in 
Fig.  11  illustrates  that  the  correction  term  is  bounded 
everywhere.  In  fact  the  spikes  are  completely  absent  at 
the  locations  of  maximum  and  minimum  amplitudes 
(and  zero  strain  rate);  at  these  locations  the  correction 
L\  is  identically  zero. 

From  the  analysis  and  computations  of  this  section 
we  conclude  that  the  power  law.  Cross,  Carreau, 
Prandtl-Eyring,  and  Powell-Eyring  forms  are  problem¬ 
atic  as  characterizations  of  a  strain-rate  thinning  fluid  in 
flows  experiencing  locations  or  times  of  low  strain  rates, 
since  they  cannot  model  observed  slender  oscillating  jet 
behavior:  The  first  order  corrections  to  the  leading  order 
solution  are  either  unbounded  (for  the  power  law  and 
Cross  forms)  or  large  (for  the  Carreau,  Prandtl-Eyring, 
and  Powell-Eyring  forms)  at  the  locations  where  the 
strain  rate  is  passing  through  zero;  these  isolated  spikes 
indicate  that  if  the  fluid  viscosity  was  in  fact  governed  by 
the  assumed  form,  an  oscillating  jet  could  not  be 
established,  in  opposition  to  physical  reality.  Our 
proposed  Zhou  form,  on  the  contrary,  yields  uniformly 
small  first  corrections,  and  hence  passes  the  test  of  being 
able  to  model  oscillating  jet  behavior. 


Conclusions 

Oscillating  jets  with  an  elliptical  cross-section  provide 
accurate  means  for  material  characterization  of  other¬ 


wise  difficult  properties,  namely  dynamic  surface  ten¬ 
sion  and  elongational  viscosity.  Here  we  have  focused 
on  the  relevance  of  oscillating  jet  models  that  accom¬ 
pany  this  experimental  technique  for  the  specific  case  of 
strain-rate  thinning  and  strain-rate  thickening  viscous 
fluids.  We  have  shown  analytically  and  illustrated 
numerically  that  standard  viscosity  forms  (the  power 
law.  Cross,  Carreau,  Prandtl-Eyring,  and  Powell- 
Eyring  forms)  for  strain-rate  thinning  inelastic  fluids 
break  down  when  applied  to  oscillating  jet  flows  in 
which  gravity  is  a  weak  effect,  precisely  because  of 
behavior  of  these  forms  near  zero  strain  rate.  We  then 
propose  a  new  form  that  can  fit  data  at  least  as  well  as 
the  Cross,  Carreau,  Prandtl-Eyring,  and  Powell-Eyring 
forms,  and  is  in  addition  consistent  with  the  asymptotic 
model.  These  results  offer  a  new  model  for  material 
characterization  experiments  on  strain-rate  thinning 
fluids. 
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